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Abstract 

This paper treats the global existence question for a collection of 
general relativistic collisionless particles, all having the same mass. 
The spacetimes considered are globally hyperbolic, with Cauchy sur- 
face a three-torus. Furthermore, the spacetimes considered are iso- 
metrically invariant under a two-dimensional group action, the orbits 
of which are spacelike two-tori. It is known from previous work that 
the area of the group orbits serves as a global time coordinate. In the 
present work it is shown that the area takes on all positive values in 
the maximal Cauchy development. 

1 Introduction 



Global existence results have been obtained [TJ El HI El EH Ej for various 
classes of T 2 symmetric spacetimes - solutions to Einstein's equations with 
a closed (compact without boundary) Cauchy surface isometrically invariant 
under an effective group action, U(l) x U(l). A key feature of these results 
is that the area of the two-dimensional group orbits serves as a global time 
coordinate. In [TJ El Ej the global existence result is obtained directly in 
terms of a coordinate system with areal time coordinate. In these cases the 
"areal coordinates" are also "conformal coordinates" and the area tending to 



1 



zero signals the end of the maximal Cauchy development. 1 In the other cases 
|U 03 [7j , areal coordinates and conformal coordinates are distinct, and the 
global existence argument in the contracting direction is not made directly 
in terms of areal coordinates, but rather in terms of conformal coordinates. 
(The quotient of the spacetime by the symmetry group is conformally flat. 
In conformal coordinates the conformal factor is given explicitly.) The re- 
sults in the second group leave open the possibility that any conformal time 
coordinate is unbounded in the contracting direction and that the area tends 
to some positive number (rather than all the way to zero) at the end of the 
maximal Cauchy development. 

In |B], vacuum T 2 symmetric spacetimes such that the twist quantities 
associated with the two commuting spatial Killing vector fields do not both 
vanish are considered. Areal coordinates are used directly in the contracting 
direction to show that the area of the T 2 symmetry group orbits tends to zero 
at the end of the maximal Cauchy development (as long as the spacetime is 
not flat - one can choose two Killing vector fields tangent to a Cauchy surface 
in a spatially compactified flat Kasner spacetime such that the twist quan- 
tities are not both vanishing and such that the area of the T 2 group orbits 
tends to a positive number at the end of the maximal Cauchy development). 
The argument used in [H] to thus sharpen the vacuum global existence result 
previously obtained in [3] carries over straightforwardly to sharpen in the 
same way the global existence result for spacetimes with Gowdy symmetry 
(i.e., vanishing twist quantities) and a nonvanishing magnetic field orthogo- 
nal to the group orbits obtained in appendix C of 0. In the present work, 



the general strategy of jH], using areal coordinates directly in the contract- 
ing direction rather than conformal coordinates, is used to so sharpen global 
existence results previously obtained in [3J [7] for Einstein- Vlasov initial data 
on T 3 with T 2 symmetry. 

The Einstein- Vlasov system models a general relativistic collection of col- 
lisionless particles. For an introduction to the Einstein- Vlasov system, see 
[Jill EEH U2] • In the present work, as in jU [7] , massive particles are con- 
sidered, all having the same mass. In section |21 the areal coordinate metric 
is written out, a global existence result from [7] is recalled, and the Vlasov 
matter distribution function is introduced. Some preliminary results per- 

1 The spatial topology is the three- torus in ^0 E3 IZ| and a global time coordinate 
can be set equal to the area of the group orbits. In 0, other spatial topologies are 
considered, and the relationship of the area of the group orbits to the time coordinate is 
specified. 
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taining to the Vlasov matter, which are key to the presence of Vlasov matter 
insuring that the area of the T 2 symmetry orbits goes to zero, are obtained. 
In section El Einstein's equations are written out. In section 0J bounds on 
the areal metric functions are obtained in the contracting direction which are 
sufficient to allow application to the contracting direction of the techniques 
used in j2HU[7] for extending the spacetime in the expanding direction. These 
techniques involve integrating quantities related to the metric functions along 
the light cones of the quotient spacetime and integrating quantities related 
to the Vlasov matter along timelike geodesies. The main theorem appears in 
section 

There does exist an exceptional possibility for T 2 symmetric spacetimes, 
characterized by the spacetime gradient of the area of the symmetry group 
orbits vanishing on the closed Cauchy surface, which implies that the area of 
the orbits is everywhere constant, and so cannot serve as a time coordinate. 
In this exceptional case, the spacetime is flat, either T 3 xR Minkowski space- 
time or flat Kasner. The presence of a nonvanishing twist quantity rules out 
this exceptional possibility [2 , as does the presence of Vlasov matter (see 
the appendix). 

2 Preliminaries 

One of the results in [7j is that (except for the special vacuum solutions with 
symmetry group orbits of constant area mentioned in the introduction) the 
maximal Cauchy development of T 2 symmetric C°° initial data on T 3 for the 
Einstein- Vlasov equations 2 is isometric to (to, oo) x T 3 with C°° metric, 



for some nonnegative number to- All the metric functions are periodic in 9 
with period one and independent of x and y so that X = d x and Y — d y are 
commuting Killing vector fields. The area of the symmetry orbits is t. (The 
result in 0] is similar, but with an additional assumption made, that G and 
H vanish.) 

The goal of the present work is to show that to = 0, that the area of the 
group orbits goes to zero at the end of the maximal Cauchy development, 

2 The area of the group orbits need not be constant on the initial data surface. The 
hypotheses concerning the Vlasov matter are stated later in this section. 



9 = 



e 2iu - u \-adt 2 + d6 2 ) + e 2U [dx + Ady+{G + AH) d0} 2 




(1) 
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in the contracting direction. The strategy for showing that to = will be 
to show that, given a C°° solution to the Einstein- Vlasov equations with 
metric (JIJ on (to, oo) x T 3 , and a nonvanishing number of Vlasov particles 
(see below), then for any t p G [to, oo) such that t p > 0, the spacetime extends 
to t p . It will be convenient to fix t p G [to, oo) such that t p > 0, and to fix 
some arbitrary tj G (t p , oo). 

There are two twist quantities associated with the commuting Killing 
vector fields, J = e w<JT X^Y p V a X T and K = e WCJT X p Y p V a Y T . If the orien- 
tation is such that {d t , dg, d x , d y } is right handed, then these twist quantities 
are related to the metric functions as follows, 3 

f -2u+4U +Z e ~2v 

J = T=^(G t + AH t ) and K = AJ —H t . (2) 

v« v« 

Consider massive Vlasov particles all of the same mass. Without loss of 
generality, the mass can then be normalized to one. Let the velocity vector of 
the Vlasov matter be future directed (with the future designated as increasing 
t). Let the quantities, v^, denote the components of the velocity vector in 
the frame 

{dt, d6, dx + Gd6, dy + Hd6}. (3) 
With this choice of frame, i>2 and v 3 are constant along geodesies, and from 

v = -^ae^~ 2U + avj + ae 2v ~ w v\ + at- 2 e 2 »{v 3 - Av 2 ) 2 . (4) 

The Vlasov matter distribution function, /, is in general a nonnegative func- 
tion of position and velocity (or momentum). Suppose that C°° initial data is 
given for the Einstein- Vlasov equations 4 on some T 2 symmetric initial three- 
torus, with a nonvanishing number of Vlasov particles (see the paragraph 
following lemma 121 below) , with the matter distribution function invariant 
under the action of the symmetry group and with the support of / bounded 
initially. From [7j it is known that / is C°° on (to, 00) x T 3 x R 3 , and the sup- 
port of / is bounded at each t G (to, 00). The Vlasov equation is that / should 
be constant along timelike geodesies. In coordinates [t,6,x,y,vi,V2,v^), the 
Vlasov equation takes the form 

3 With this choice of orientation, K as defined here is related to K as defined in |3] and 
[Hj by change of sign. J = in and [HJ- 

4 See ^3 Ell d] for statements of the Einstein- Vlasov initial value problem. 
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df_ dv df f dv 
dt ~ dvi d6 I 86 

\/™?" its atm as Vae 2u ~ 4U T \df 
+ }L ^( K - A J )( v 3 ~Av 2 ) + Jv2j-^, (5) 

with f x = 0, f y = and / periodic in 9 with period one. 

Lemma 1 There exists a number G\ > such that f < C\ on (t p , ti] x T 3 x 
R 3 . 

Proof: This follows from / being constant on timelike geodesies. □ 
The second lemma follows from the continuity equation, V M iV M = 0, 
which holds generally for Vlasov matter, with the matter current density. 
In the frame 0, the matter current density is 



N » = f ^Ld Vl dv 2 dv 3 . 

t JR? \Vn\ 



la, f fv^ 1 

IB? \vq\ 

Lemma 2 The quantity C 2 = Is^Ur? f dv\ dv 2 dv$) d9 is constant in time. 

Proof: This lemma can be obtained directly from the continuity equation. 
Similarly, it can be obtained from equation (JHJ), as follows. 

dC 2 r ( r f dv df rdv ^fae 1 " 



dt 



r fr (dv df rdv yae" 



-Jn 



„ — ? dv\ dv 2 dv3 I d9, (6) 
t J dv 1 J / 

= 0. (8) 

Integration by parts leads from line (JHJ) to line (jZJ. □ 
The property that the number of Vlasov particles is nonvanishing is equiv- 
alent to C 2 0, which is assumed throughout the rest of this work. 

Since / doesn't grow arbitrarily large, since the integral of / is a constant 
in time, and since /, v 2 and v 3 are all constant along timelike geodesies, 
f s i {Jr3 f\vi | dvidv 2 dv 3 )d9 does not become arbitrarily small on (t p , ij, as 
shown by the following lemma. 
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Lemma 3 There exists C 3 > such that J s i (J R3 f\v i| dvidv 2 dv 3 )d9 > C 3 
on (t p ,ti\. 

Proof: Fix positive numbers v 2 and v 3 such that the support of / at ti is 
contained in S 1 x R l x [— ^2,^2] X [—^3,^3]- At any t G (t p ,U], the support 
of / is contained within this same set. This is because /, v 2 and v 3 are 
constant along timelike geodesies, so that if / vanishes for all 9 and v\ at 
some (ti, v 2 , v 3 ) then / vanishes for all 9 and v% at (t, v 2 , v 3 ) for all t G (t p , tj\. 
Fix 5 > small enough so that b = C 2 — 85v 2 v 3 Ci is positive. From the 
preceding paragraph and lemma [H 

L ( L ( 1-6 f dvi ) dv2 dv3 ) de - m ^ Ci - (9) 

Now using lemma |2l 

/ ( f I / / <foi J d^d^ > b. 

JS 1 \JR 2 yJ^yS J J 

To obtain lemma El calculate 

J (^J^f\v\\dvidv 2 dv 3 jd9 = J (^J 2 (^J f\v\\ dv-^j dv 2 dv 3 S j d9 

+ LU# {L>s flvildvi ) dv * dv *) de 

> 5 (^J r2 n | JdvA dv 2 dv^j d9 

> 5b, (10) 

and let C 3 = 5b. □ 
In the course of the proof of the main result, a coordinate transformation, 
(t,9,vi,v 2 ,v 3 ) — > (s,4>,W\,w 2 ,w 3 ), will be used, with t = s (letting s, = U 
and s p = t p ), v 2 = w 2 , v 3 = w 3 , 89/ ds = —dv^/dvi, 

dv, dv ^ 2u \ (K-AJ)(v 3 -Av 2 ) e~ iU Jv 2 \ 

97 = w + ^ e 1 * + ^^J' 

9(si, (p,Wi,w 2 ,w 3 ) = (ft and t>i(sj, <fi, W\, w 2 , w 3 ) = u>i. The new coordinates, 
(0, wi, w 2 , w 3 ), are constant on timelike geodesies and therefore f s = 0, which 
is why the coordinate transformation is useful. 
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Lemma 4 The Jacobian determinant of the transformation (t, 9, v%, v 2, v 3) 
(s, (f), w%, w 2 , w 3 ), equals one on (s p , Sj] xS'x R 3 . 



Proof: The Jacobian determinant of the coordinate transformation is 

89 dvi 89 dvi 
^ dcj) dw\ dui\ d<p 

On {si} x S 1 x i? 3 , 77 = 1, and on (s p , Sj] x S 1 x i? 3 , 



d_/de\de_dvi_ 8 fdO\ dvi dvi 
89\dsJ d<p dw\ dv\ \8sJ dcj) dw\ 
89 8 90 8 /8v 1 \8v 1 

8<j> 89 \ 8s J dw\ 8<j> 8v\ \8s J 8w\ 
8 (d0\ 89 dv 1 __ 8 (d6\ 8v x dv 1 
89 \ds J dwi 8<j> dvi yds J dwi 8cj> 

_d9_d_fdvi\d9_ _ ^^/8vi\8vi 
dw\ 89 \ 8s J dcj) 8w\ 8v\ \ 8s J 8<j> 



From 



8_f89\ d /dv x \ _ 
89\8s) dvi V ds J 

and cancellation of terms in ([lip, it follows that r) s — 0. □ 
In general, coordinates that are constant on timelike geodesies are valid 
on the maximal Cauchy development of an initial data surface. Here one 
can see this also directly from lemma |U combined with the smoothness of the 
metric functions appearing in on the maximal Cauchy development. 

3 Einstein's equations 

With choice of frame (JHJ), the stress-energy tensor due to the Vlasov matter 



is 



« f fv^Vp 



t Jrz \v \ 
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Einstein's equations, G w = T MP , for the metric (JTJ) with source (f*HZ|) are 
satisfied if 

n = tlu 2 + aU 2 + ^(A 2 + aA 2 )\ + ^^J 2 



ae 2 » 



|/3 (K-AJ) 2 + yfaj f\v \ d Vl dv 2 dv 3 , (13) 



2t (u t U e + ^pA t A^j - ^ - J r3 f v 1 dv x dv 2 dv 3l (14) 



at ae 2v ~ AU J 2 ae 2v (K - A J) 2 

a ~ t t 3 

_ 2a3/ ^ I Ke-™ + e-«vl + t-^-Av 2 ?) ivi 
Jr3 I v | 

J t = —2a I —. 1 2 dv i (if 2 (if 3, (16) 



B 3 | Vo| 

J e = 2 fv 2 dv 1 dv 2 dv 3 , (17) 

JR 3 

i^t = —2a / ^ 1 3 dv\ dv 2 dv 3 , (18) 
Jrz I v | 

K e = 2 fv 3 dv 1 dv 2 dv 3 , (19) 

J ij 3 

rr rr ^ , ^ , ao£/fl , e *" ( a 2 /i 
U tt = oJJqq ~ A t — aA b 

t 2a 2 2t 2 V 
ae 2 - 4!/ t2 a 3 / 2 e 2 ^ /■ /(l + 2e~ 2 % 2 h , , 

^ = a ^ + ^ + ^ + ^-4(^C/ t -a^) 

2i/-4tf 

+ - J(X-AJ) 

+ 2a3/V ^ / ^-^d^du,^. (21) 
t ,/fi 3 |v | 

Since J and K are periodic in equations (fTTj) and (fE^) impose requirements 
on /. Namely, J s i J R3 fv 2 dv\ dv 2 dv 3 d6 = and J s i J R3 fv 3 dvi dv 2 dv 3 d9 = 
0. Similarly, equation (fT4^) imposes the condition 



4V 



2t ( W« + ^A t A^\ -y/a~J R Jv x dv x dv 2 dv 3 j d# = 0. (22) 



All three of these conditions are preserved by the evolution. That this is 
so is equivalent to the integrability conditions for equations (fTEj) and (fT7|). 
equations (fTHj) and (fTUj) . and equations (fTB^l and (HU), respectively. 

It is useful to also have at hand an equation which can be derived from 
(0) and (fTB")) - (}2~T]) or, alternatively, can be obtained directly from Einstein's 
equations, G„ p = T MP , 



4 Extending the spacetime 

Next a number of lemmas are presented which give sufficient control of the 
metric functions in the contracting direction so that, later in this section, a 
slight modification of the combination light cone/timelike geodesic arguments 
used in j3] [7j for the expanding direction can be applied to the contracting 
direction. These arguments give control of U and A and their first derivatives 
and also the support of / (proposition QJ, followed by control of the first 
derivatives of / (proposition 12}. C°° bounds on all the metric functions and 
on / follow from light cone/timelike geodesic arguments and the Einstein- 
Vlasov equations. 

Lemma 5 For any number b, ^fae 2u+bU is bounded on (t p , ti\ x S 1 . 



u tt = av de + 




(23) 



Proof: 





(24) 



> 0. 



(25) 
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So on (t p , ti] x S 1 , 



b 2 /8 



2v{ti,e)+bU(Ufi) 



62 /8 

□ 

Lemma 6 For any positive number r and any number \, 

a r/2 e 2ru + XU A 2 

is bounded on (t p ,ti\ x S 1 . 

Proof: An SL(2, R) transformation of the Killing vectors relates the quantity 
under consideration here to quantities of the form considered in lemma |5J 
Consider tilded coordinates, defined by = d t , dg = do, 

X = d i = aX + bY, 

and 

Y = 8y = cX + dY, 

with a, b, c and d constants such that ad — be = 1. Then the form of the 
metric is unchanged from (£[]), from which it follows that the tilded metric 
functions satisfy "tilded" equations (f*T3|) - (j2HJ). The relation between the 
tilded metric functions and the untilded metric functions is 

e 2ij = e 2U (a + Ab) 2 + t 2 e~ 2U b 2 (26) 

e 2 °A = e 2U (a + Ab)(c + Ad)+t 2 e- 2U bd (27) 

a = a (28) 

ae 2 *- 20 = ae 2v ~ 2U (29) 

Note that t = t and 9 = 9 and now consider a = and b = 1. Choose a 
positive number q < r. From equations (|2T)j). (|2~%j) and (|2T?j) . 

^q/2 e 2qu+2(l-q)U _ Q , 9 /2 e 2 ? ^+2(l- 9 )£/^2 + Q ,g/2^2 e 2 g! /-2(l+ (? )C/ _ ^gg^ 

Since the tilded metric functions satisfy the same equations as the untilded 
metric functions, the left hand side of equation (j30j) is bounded on (t p , ti] x S 1 
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from lemma 03 Since the second term on the right hand side of equa- 
tion PU) is positive, a i/^ e ^+^-i)u j£ is bounded on [t^U] x S 1 . Since 
a (r-q)/2 e 2(r- q )v+(\-2+2q)u ig bounded on (t p ,ti) x S 1 from lemma lemma© 
follows, by taking the product of these two bounded quantities. □ 
Next the monotonic quantity from [7] will be recalled. In [7| this quantity 
was called E(t). The notation here follows [H]. Let 

set) = ( ^de. (31) 

Lemma 7 dS/dt < 0. 
Proo/: 

1 7 hfe + TaV^ + V ^ J 2 + ^~— (i^-AJ) 2 



dt ^ U Wo 4t JY y / 2t 3 t 5 



Equation ()32|) follows from the evolution equations (|15|) and ()23jl and from 
the vanishing of the integral over the circle of the derivative with respect to 9 
of any C 1 function on the circle. A nonvanishing number of Vlasov particles 
(C2 7^ 0) insures that dS/dt is strictly negative. □ 

Lemma 8 S(t) is bounded on (t p ,ti\. There exists a number, £ p , satisfying 

S v = lim£(£). 

tit p 

Proof: From equation (J32*Jl . dS/dt > —A£/t. So for any tk G (t p ,ti\, 

S(t k ) <£{U)+ P^P-dt. (33) 
Jt h t 

Applying Gronwall's lemma to this inequality (as suggested on p 353 of |lj), 
we obtain 

( 34 ) 



for any tk € (t p , tj. Thus S(tk) < S(ti){U/t P Y on (t p , tj). This bound, 
together with the monotonicity of S(t), guarantees that S p , as defined in the 
statement of the lemma, exists. □ 
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Let 



P = v + 



In a 



(35) 



2 



The substitution of e 13 for \fae v will be made freely throughout the rest of 
this section. This is a convenience, because of the order in which bounds 
are obtained. The advantage of j3 over v is that ag does not appear in the 
expression for f3e obtained from equation (jHJ), 



Lemma 9 There exists a positive number, C±, such that j s i \j3g\ d6, j s i \Ug\ d6, 
J s i e 2U \A g \ d6, J sl \Je\d6 and J s i \K e \d6 are each less than C± on (t p , tj\ . Thus, 
the difference between the spatial maximum and minimum values of each of 
the quantities, (3, U, J and K on S 1 (at any fixed time t) is less than CV 

Proof: Using that a and S are increasing with decreasing t, 




(36) 




(37) 



(38) 



(39) 



(40) 



(41) 



on (t p ,ti). 



□ 



Lemma 10 min5 1 a(t,^) is bounded on (t p ,ti\. 



Proof: For any t G (t p ,ti], (see equations ()31|) and f(13jl ) 




(42) 



Now considering equation (jlj), 
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min a(t, 9) 



Therefore, 



Si 



From lemma El 



with C, > 0. So 



From lemmas [7| and |Hl 



S 1 JR? 



f\v\\ dvi dv 2 dv 3 d6 < t£{t). 



mina(t,#)C 3 < t£{t), 



t£(t) 



min a(t, 9) < ( , 



on (t p ,ti). 

Lemma 11 There exists 9 G Si such that a(t,9) is bounded on (t p ,ti] 



(43) 
□ 



Proof: Suppose there existed no 9 such that a(t, 9) were bounded on 
(t p , ti]. Since a(t, 9) is increasing with decreasing t, this would contradict the 
previous lemma. □ 



Corollary 1 On (t p ,ti\, 



da 



t a(a. 



is bounded. 

Proof: In a(t, 9) is bounded since a(t,9) is bounded and increasing with 
decreasing t. □ 
In the next four lemmas control of the metric functions is obtained which 
is sufficient for using the light cone/timelike geodesic arguments directly in 
terms of areal coordinates in the contracting direction. 



Lemma 12 For any number h, e^ +6C/ is bounded on (t p ,ti\ x S 1 . 
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Proof: It follows from lemma along with a being bounded on {t p ,tj\ x {9} 
that e^ +bU is bounded on (t p , ij] x {9}. (See equation fj35j) for the definition of 
(3.) That e^ 3 " 1 " 617 is bounded on (t p , ti\ x S 11 then follows from the boundedness 
of the difference between the maximum and minimum values of (3 and U on 
S 1 (see lemma EJ). □ 

Lemma 13 For any positive number r and any number b, 

e rP+bU A 

is bounded on (t p ,ti\ x S 1 . 

Proof: It follows from lemma El along with a being bounded on {t p ,tj\ x {9} 
that e r/3+bU A is bounded on (t p ,ti\ x {#}. Integrating (e r ^ +bU A)e along at = 
constant path from (t, 9) to (t, 62) leads to 



r/3+w/ A)(t,0 2 



rl3+bU 

e 2 



A)(t,9) 



+ J_ ^ +bU A{r(3 e + bU e ) 
+e rP+(b-Vu e 2U A \ ddj 

?fWA)(t,9 2 )\ < \{e rP+bU A)(t,9)\ 



(44) 



+ 



02 



e^ +bU \A\{r\(3e\ + \bU e \) 



< c 5 + 



02 



> rf5+bU \A\{r\(3e\ + \bU e \)d9 



(45) 
(46) 



on (t p ,ti\, for some number C 5 > (see lemmas IP21 and Ej) ■ From Gronwall's 
lemma, 

\(e rP+bU A)(t,e 2 )\ <C 5 exp| 



(r|ft| + |&E/*|)d0 



(47) 



which (using lemma EJ) determines a bound for e r/3+6C7 A on tj\ x 5 1 . □ 
Lemma 14 

max^^-^^J^, £)) 2 ] da 

is bounded on (t p ,ti]. 
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Proof: From lemma El 



max[e 2 ^- 4 ^J(M)) 2 ] da 



< { * e 6C *e 2 ^- 4U ^\j(o-, Of + 2C,\J(a, 6)\ + Q 2 )da. 



(48) 



The right hand side of inequality (J48j) is bounded from corollary Q the evo- 
lution equation (|15|) and lemma IT21 □ 

Lemma 15 

/ max [ e WW)(K(a, 6) - A(a, 6) J (a, 6)) 2 ]da 

J t $ 

is bounded on (t v ,tj\. 

Proof: Integrating [ePiK — AJ))g along a t — constant path from (t,9) to 
(t, 9 2 ) leads to 



{e?\K-AJ\)(t,0 2 ) < (e?\K-AJ\)(t,i 



02 



J\K-AJ\\p e \ 



/ f\v 2 \dvidv 2 dv 3 

JR. 3 

+e^A / f\vs\dvidv 2 dv3 

JR? 



+eP- 2U \J\e 2U \A \jd6 

< (eP\K - AJ\)(t,6)+C e 
+C 4 max(e /3 - 2C/ |J|) 

5 1 



(49) 



+ 



e 2 



e p \K-AJ\\p e \d0 



(50) 



for any (t,9 2 ) <G (t p ,ti) x S 1 , with 



C e = C 2 (v 2 sup (e?)+v 3 sup (ePA)). 

(tp ,t{\ X S (tp X S 



From Gronwall's lemma, 

(eP\K-AJ\)(t,9 2 ) < e 



Is 2 \M dd 



{(eP\K-AJ\)(t,d)+C 6 



+C 4 max(e /3 - 2C/ |J|)|. 



(51) 
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Using lemma El 

mwc[(e?\K-AJ\)(t,6)] < e° A Ue^\K -AJ\)(t,0) + C 6 



+C 4 max(e /3 - 2C/ |J|)|. (52) 



Therefore, 



max 

3 2c 4 / - j e m°fi)( K ^ o) - A(a, 9) J(a, 9)) 2 
+2C 6 (e^\K-AJ\)(a,9] 



+2(e l3 \K - A J\)(a, 9)C 4 max{e^ 2U \ J|) 

s 

2 



C 6 + C 4 max(e /3 - 2l/ |J|)) daj. (53) 



Using that 



2(e l3 \K - A J\)(a, 9)C 4 max(e /3 - 2C/ | J\) 

< {e 2 \K - A J) 2 )((T, 9) + C\ max(e 2/3_4C7 J 2 ), 

the right hand side of inequality (J55j) is bounded from corollary Q the evo- 
lution equation (JTHJ) and lemma ITU □ 

Proposition 1 The junctions U and A and their first derivatives are bounded 
on (t p ,ti\ x S 1 . The support of f is bounded on (t p ,ti\. 

Proof: The proof of this proposition (which has to do with the contracting 
direction) is a slight modification of the argument used for the expanding 
direction in [3J and [7J. Let 

1 e iU 
E = -{U 2 + aU 2 + —(A 2 + aA 2 e )}, (54) 

P = ^(U t U e + —A t A e ). (55) 

Since E ± P are both the sum of squares, \P\ < E. Set 

(d t TV^d e )(E±P)=L ± . (56) 
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A straightforward calculation shows that 



+ (U t ±^U e )\-^J 2 

+— (A ± ^)|^- J(# - A J) 

M^^^U. (57 ) 



Let u\ = \foLV\. 



a v I 



e 2/3 



+ f ?f>-W(p e _ 2 ^ 2 + _( V3 - Av 2 )[(v 3 - Av 2 )(3e 
-A dV2] \ + 2 ^ ^ K - AJ ^- A ^ + e^J V2 ) (58) 



The advantage of u\ over v\ is that does not appear in equation (f58|) . 
Let Mi(t) = suplv^ajvil : f(t,9,vi,v 2 ,vs) ^ for some (t, t> 2 , t> 3 ) G [t, tj] x 
S 1 x i? 2 }. Analogously to the expanding direction in [3] and [7j (but here 
considering the contracting direction), the strategy to prove the proposition 
is to show that max s i E(t, 9) + u\{t) is bounded on (t p , tj\. Let 

\r,u\ - S max s 1 °) + «?(*) if max s i S(t, 9) + w 2 (t) > 1 

nrj ~\l ifmax 5 iE(t^)+M 2 (t) < 1. l ° yj 

If 

*(t) <C 7 + J* h(o)V (a) In* (a) da (60) 

on (t p ,tj], for some number C 7 > 1 and function /i(t) > 0, it follows (since 
* > 1) that 

h(t)*(t)to*(t) 

(C 7 + J* fc(<r)tf (<r) ln*(a) dcr) ln(C 7 + J** h{a)^{a) ln^(cr) da) ~ 1 J ' 
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Integrating both sides of this inequality and then exponentiating twice and 
using inequality (jBDJl one finds that 

exp ( f' 1 h{a) do) 

*{t) < C 7 V 1 (61) 

on (t p ,ti\. It will be shown that \l/ satisfies inequality (jBTIj) with the quantity 
J t * h{a) da bounded on (t p ,ti), so that inequality (JBTj) determines a bound 
for \l/ on (t p , ti]. 

First consider maxgi i?(t, 6 1 ). Consider an arbitrary point (tk, 8k) € (t p , ti] x 
S 1 . Let 7 ± be the integral curves of d t ± y/ade starting from the point (t&, 6^) 
and extending to the surface t = ti. Let (ti,9±) denote the endpoints of j± 
lying on the surface t = ti. Then 

/ (d t - v ^d )(E + P)+ [ (d t + yfed e )(E-P) = I L + + f L_. (62) 

J'Y— ^7+ *^7— *^7+ 



So 



L — L 



+ h 

7+ 



< E{t u 6 + )+E{t h e-) 

All \ L -\+f \ L + \\- (63) 



2 



7+ 



Consider L± (see equation (|37jl). The terms in the brackets satisfy 
e 2 ^ 4f/ T2 v^e 2 ^ /■ /(l + 2e" 2 % 2 ) , , , , a t , 

e2 "' 2U J(K-AJ ) + ^L[ i^^ dvM <\^\. (65) 



2t 3 v ' t 2 Jrz \v \ '2ot 

The sum of the factors in front of the brackets satisfies 

p 2U AU i 

C/ t ±VSC/ + — (A t ±^A e ) < (U t ±^U e ) 2 + —(A t ±^A e ) 2 + - 

< 2{E±P) + l - 

< ±E+ 1 -. (66) 
18 



Thus 

\L ± \<\^(2E +T + -) + T (67) 

The first two terms in a t /a (see equation (|15|l ) are treated in lemmas PUP 
and El Changing variables from to u±, the third term satisfies 

Jrz \v \ 
< 2e^ / /(l + e-^ + r^fa-^)-) 



t 2 y-fi! + e ~ 2 P +2U ui 

+£ 72-(^ + \A\v 2 ) 2 )v 2 v,(e~ l + e^lnfa + ^ e 2 ^ + u 2 )) . (69) 



One of the factors, e f3 ~ u , appearing in front of the integral in line (J68|) is 
used to control the terms appearing in the numerator and the other has been 
absorbed into the integral in the following line, to control it, since it has not 
yet been shown that e@~ u is bounded away from zero. Now using lemmas 
and ITSl and inequalities (ffil3|). (joTj) and (pj9*|) . 



maxE(t k ,6) < 2m&xE(U,6) 

+ 



51 



tk 



I t 

max s i [e 2 ^%K(t, 6) - A(t, 6) J(t, 9)) 2 } 
+ T 3 

+C 8 (e" 1 + eP~ u mtmax^x, e^ u })\{2E + ^ + i-) 

+^]dt, (70) 

for some positive number Cg. 

Next consider -u 2 . Integrating along geodesies using equation (jSHJ), treat- 
ing «t/a as above, and using that \ui\ < \vq\, that 

Va\/3e\< 2tE + Au\v 2 v z C u (71) 
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that y/a\U e \ <2E + 1 and that y/ae 2U \A e \/t < 2E + 1, 

max 5 i [e 2/3 (*- e ) 0) - A(t, 6) J(t, 6)) 2 } 



+ 



t 3 



+C 8 (e" 1 + e 73 " 17 ln(max{fti, e^"^}) | (n x (t)) 2 

+C 9 (£ + ft? + 1)| maxe^- 2U + (maxe^ 417 )^ 

[(max s i e /3 )v 3 + (max s i e^ADf^] 2 
H i 2 

| [(maxgi e 2 ?- 2U )v 3 + (max 5 i e 2/j - 2C/ |A|)^ 2 ]t; 2 
max^e^l-fT — A J|)[(max5i(e 73 )f;3 + maxsi(e^|A|)f;2] 



f 3 

(max 5 i e 73 " 217 ) max 5 i (e 73 " 217 1 J\)v 2 



))dt (72) 



for some positive number Cg. From ()70|) and ()72|). and using lemmas ITU and 
1131 it follows that ty(t) satisfies inequality (|50|l with /i(cr) da bounded on 
so that inequality (pTTj) determines a bound for \1/ on (t p ,ti]. That 
the support of / and U t are bounded on (t p , ti] x S 1 follows immediately. 
Therefore U is bounded and from this follows the bound on A t and therefore 
A. The bounds on the derivatives with respect to 6 follow from the bound 
on E(t, 9) and the fact that a is increasing with decreasing t. □ 

Corollary 2 In a is bounded on (t p ,ti] x S . 

Corollary 3 J, K and their first derivatives are bounded on (t p ,ti] x S 1 . 

Corollary 4 The functions a t , (5, (3 t , (3q, v and v t are bounded on (t p , t^xS 1 . 

Proposition 2 The first derivatives of f are bounded on (t p , ti] x S 1 x R 3 
and ctg is bounded on (t p , ti] x S 1 . 

Proof: The appropriate quantities have now been bounded so as to allow 
this argument to proceed as for the expanding direction in 0] and [7]. The 
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orthonormal frame used in 0] and [7j is such that the components of the 
velocity vector are 

v° = -e-P +u v , (73) 



v 



i = ^e^ +u v 1 , (74) 
v 2 = e - u v 2 , (75) 

v 3 = ^(v 3 -Av 2 ). (76) 



t 



Note that v° = y 1 + (-0 1 ) 2 + (v 2 ) 2 + (v 3 ) 2 and that the support of / is 
bounded in terms of these variables also. Let Q = (s,(f),Wi,W2,Ws), R = 
(t,9,v 1 ,v 2 ,v 3 ) and S = (t, 9, v 1 , v 2 , v 3 ). Let Q° = s, Q 1 = 0, Q j+t = Wi and 
analogously for coordinates R and S. Let the indices J, K and L take on 
values from to 4. The goal is to determine bounds for df/dR J . 
The first step will be to determine bounds for 

df df dQ K 
dS 7 ~ dQ K dS J ' 

Since df /dQ K is constant along timelike geodesies, for the first step it is 
enough to determine bounds for dQ K / dS L . Since the determinants of OS/ dR 
and dR/dQ (see lemma |3J) are bounded away from zero, the determinant of 
dS/dQ is bounded away from zero, so it is enough to determine bounds for 
dS L /dQ J . Let the indices a, 6, c and d take on values from 1 to 4. Since 
8S L /dQ° (see equations (JSTj) - (JH3J) below) and dS°/dQ J are bounded, it is 
enough to determine bounds for dS a /dQ b . Let 



F 2 

a 




(78) 

+ -. (79) 



^.3 e 
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(F b a is the same as (\I/, Z l ) in jU E].) It is straightforward to see that F£ = 
U®dS c /dQ b , with II" and its inverse bounded. Calculation of dF^/ds, using 

t s = 1, (81) 
«. = ^ <«) 

-(A - c/,)V5«° - (A - £/, - ^)«' - v^c/» ( '" 3)2 r„ ( " 2)2 

— -^--(i^-Ajy, (83) 



v 2 = "-^ 



0? 



t " t)° t t 3 
f/ t {) 2 - v 7 ^^-, (84) 

(C/ t - V + v^C/^ - V(^ 2 + V^^). (85) 



shows that dF^/ds = K a c dS c /dQ b , with bounded. (All terms appearing 
during the course of the calculation for which bounds have not yet been deter- 
mined are canceled by other terms, with use of the evolution equations ([20)). 
(J2H and p5]l.) Therefore dFg/ds = A^n -1 )^ is sufficient to determine 
bounds for F£, and therefore dS a /dQ b . 

The second step is to determine a bound for a$. From the evolution 
equation (fTBTh 

ae 2 ?- 4U J 2 ae^(K-AJ) 2 
a t = — 



t t 3 
-2ate^- / /(1 + ( ' 2) n 2 + ( ' 3)2) ^^^. (86) 

Taking the derivative with respect to 9 of both sides of this equation and 
then integrating in time leads to an inequality of the form 

\oce\<C 10 + j \a e \pda, (87) 

with Cio > constant and p{t) nonnegative and bounded on (t p ,ti], so a 
bound for ag is determined by Gronwall's lemma. 
Finally, 

dJL _df_dS^_ 

dR J 8S K dR J ' K ' 

determines bounds for df/dR J . □ 
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Corollary 5 Ug, a a and a t e are bounded on (t p , ti] x S 1 . 

Lemma 16 The second derivatives of U and A are bounded on (t p ,ti\ x S . 

Proof: Let 

If e 4U 1 

E 1 = - |(t/«) 2 + a{U te f + -^U tt f + a(A e ) 2 ]| , (89) 

P 1 = V^iU tt U t0 + ^A tt A t9 \. (90) 

A bound for maxsi E l (t,6) for t G (t P ,ti] can be obtained using; Gronwall's 
lemma by again considering an arbitrary point (tk, 6k) G (t p , U] x S 1 , integrat- 
ing (d t =F y/oidg)(E 1 ± P 1 ) along 7 T (as defined in the proof of proposition |TJ) 
and considering maxima on S 1 to obtain 



max E\t k , 9) < C u + maxE 1 ^, 6) £dt, (91) 

S J tip s 



for some positive number Cn and some function which is bounded and 
nonnegative on (t p , The bound for max s i -E 1 (t, 6>) on (t p , ij, determines 
bounds for £7^, U t e, A tt and Ajg on (t p ,ti] x 5 1 . The bounds for C/^ and 
then follow from equations (J2*U|) and (|2ip. □ 

Corollary 6 TTie second derivatives of J, K and (3 are bounded on (t p ,ti] x 
5 1 . and z/ 4i are bounded on (t p ,ti] x 5 1 . 

Proposition 3 Let be any positive integer. All kth order derivatives of 
the matter distribution function and of the metric functions are bounded on 
(tp,ti\ x S 1 x R 3 and on (t p ,ti\ x S 1 , respectively. 

Proof: This follows from an inductive argument, as follows. Let n be any 
integer greater than one. Suppose that bounds have been determined for all 
derivatives of u, j3, a, U, A, J and K up to and including order n except 
for d n p/dd n and d n a/d6 n . In addition, suppose that bounds have been 
determined for dS/dQ up to and including order n— 1. Let 



d n S a 1 d"-ir a 86 39 c9 

(n) bl - bn - dQbl ; ; ; ^ + ^ 9 g 6l ■ ■ ■ dQbn , 
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with 



r 1 
r 2 



'a 

[vt ~ U t )v° 
+U t 



(93) 



v 3 ) 2 - (v 2 ) 2 







J2U 



A 



v 

v 2 v 3 



+ (U t - - ^U 6 )v x 



m 2 - (v 2 ) 2 



[v ) 2 - (v 

" 1 " 2 " 3 

V V V 



1\2 



p 3 

r 1 



\faUQV 1 



J1U 



(94) 
(95) 
(96) 



Similar to the proof of proposition |2l integrating d s F^ bi ^ bn along timelike 
geodesies determines bounds for d Q®" S QQb n , followed by a bound for d n a/d9 n 
and finally bounds for d n f/dR Jl ■ ■ ■ dR Jn . The evolution equations now de- 
termine bounds for d n v/d9 n , and the (n+ l)th order derivatives of a, except 
for <9 n+1 a/<96 |n+1 . Now bounds can be determined for the (n + l)th order 
derivatives of U and A by carrying out n different light cone arguments, with 



E n ±P n 

mr mr 




± \J~OL- 



+ 



d n - l U tt . 



d6 m dt r 

AU (Qn-X Att 



te 



d9 m dt T 



At 2 V d6 m dt r 



gn-l A 



te 



d9 m dt r 



(97) 



for each choice of nonnegative integers m and r such that m + r = n — 1. 
Finally, the evolution equations determine bounds for all (n + l)th order 
derivatives of J, K and /3, and i/ except for d n+1 u/d8 n+1 . □ 



5 Conclusion 

Theorem 1 The area of the T 2 group orbits takes on all positive values 
on the maximal Cauchy development of T 2 symmetric initial data for the 
Einstein- Vlasov equations on T 3 with a positive number of Vlasov particles 
on the initial data surface, all of the same nonzero mass. 

Proof: Suppose that t p = t . From the bound on the support of the matter 
distribution function and the C°° bounds obtained in section^ it follows that 
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there is a C°° extension of the metric functions and the matter distribution 
function to t p (satisfying the Einstein- Vlasov equations). Therefore {t p } x S 1 
is contained in the Cauchy development of any Cauchy surface in (t p , oo) x S 1 . 
This contradicts that the Cauchy development does not extend to {to} x S 1 . 
Since t p was allowed to be any positive number in [t , oo), it must be the case 
that t = 0. □ 
In [7] it is shown that there is a crushing singularity as t [ to and that there 
exists a foliation of the maximal Cauchy development by compact spatial 
hypersurfaces of constant mean curvature. 
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7 Appendix 

An argument by Chrusciel concerning vacuum T 2 symmetric spacetimes on 
T 3 (see pp 113-114 of j2]) carries over to Einstein- Vlasov T 2 symmetric space- 
times on T 3 to show that the presence of Vlasov matter guarantees that the 
spacetime gradient of the area of the T 2 symmetry orbits is everywhere non- 
vanishing. In this appendix the argument is recalled. The areal form (0) 
of the metric is only valid if the spacetime gradient of the area of the T 2 
symmetry orbits, t, is nonvanishing. However, the conformal form of the 
metric, 

g = e 2 ^- u \-dt 2 + d6 2 ) + e 2U [dx + Ady+{G + AH) Sf 

+e- 2U t 2 [dy + Hd6} 2 , (98) 

in terms of conformal coordinates (t,8,x,y), is valid even if the gradient of 
t vanishes. (Given a spacetime with metric one can set d x = d% and 
d y = dy, but to insure zero shift it may be that x ^ x and y ^ y.) The twist 
quantities are related to the conformal metric functions by 

J = -te' 2i/+AU {Gi + AHi) and K = A J - t 3 e~ 2C 'H f (99) 
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Let the frame for the velocity of the particles be 



{dt, d9, di + GdO, dy + HdO}. (100) 



Then 



v = -^e 2 »- 2U + v\ + e 2 "- iU v 2 2 + t~ 2 e 2i> {v 3 - Av 2 ) 2 . (101) 
Let h + = hf + hg and h_ = h~ t — h$ for any function, h. Let 

p AU p 2u-4U 

B = t {Ul + —A 2 + ) + - J 2 

2u r 

+ — (K -AJ) 2 + f(\v \ - v x ) dv x dv 2 dv 3 , (102) 

4:t 6 JR 3 

p iU p 2v-AU 

D = t{u * + e —Al) + e —ji 

2u f 

+ — (K -AJ) 2 + f(\v \ + v x ) dv x dv 2 dv 3 . (103) 

Note that B and D are both strictly positive since |ui| < |uo| and since, 
by assumption, the number of Vlasov particles present is nonvanishing (see 
lemmaEland the paragraph following its proof). The following equations can 
be derived from equations (5) and (6) in [7]. 

(t+) § = -B + D + t +} (104) 
(t_) s = D-u-t- (105) 

So on any t = constant surface, 

t+ (t,9) = t + (t,9 )eZ 9+{lX)dX - [* B(t,4>)e$ 0+{i ' X)dX d(t> (106) 
t4t,9) = t4t,9 )e'- f l i> - iU)dX + ( 9 D{t,<P)e-C^^ dX d<P (107) 



Thus, if t+(t,9 Q ) = for some (t,9 ), then t + (i,9) < for 9 > 9 . And 
if t-(i,6 ) = for some (t,9 ), then t_(t,0) > for 9 > 9 . This is not 
compatible with the periodicity of t± in 9, so t±(t,9o) — is not possible. 
(In vacuum and if J = and K = 0, periodicity implies that if t±(t, 9q) = 
for some (t, 9q), then t± = on the whole t = constant surface, which in turn 
implies that t± = in the maximal Cauchy development of the t = constant 
surface, and that the spacetime is flat.) 
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